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I. INTRODUCTION 

This paper extends the general results obtained in a 
recent paper by one of us, ref. [1 , for simply connected 
drums and quantum billiards to domains with a hole, 
i.e. annular drums or quantum rings. Quantum rings 
are two dimensional regions of annular shape where an 
electron is confined and possibly subject to an external 
field (magnetic or electric); as for the simply connected 
case studied in ref. [1 , only few special cases may be 
solved exactly, such as for the circular annulus, where the 
solutions may be expressed in terms of Bessel functions of 
first and second kind. However, quantum rings of general 
shape, for which exact solutions are not known, present 
interesting physical behaviors, thus justifying the effort 
of finding analytical or numerical approximations. 

For instance, it is known that the bendings of an infi- 
nite wire cause the appearance of bound states below the 
continuum threshold [3 and of localized states in quan- 
tum rings in correspondence of regions of maximum cur- 
vature [4]. In particular, Gridin and collaborators have 
formulated in Ref. [4] an asymptotic approximation to 
the modes of a quantum ring, based on the assumption 
that the ratio of the ring half- width to the radius of curva- 
ture is small. Their approach extends to ring like domain 
the classical Keller-Rubinow method O |6l 

In this paper we adopt a completely different strat- 
egy, which is both systematic and simple, and which can 
be used to obtain precise analytical and/or numerical 
approximations for the energies and eigenfunctions of a 
given quantum ring. The approach that we propose is 
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based on a generalization of the methods described in 
Ref. [Ij and allows one to solve Helmholtz equation for 
the quantum ring of arbitrary shape by conformally map- 
ping the ring to a simply connected domain. The results 
obtained in this way are very precise and prove to be use- 
ful even for the exactly solvable circular ring: in this case 
we have found an extremely precise formula for the ener- 
gies of the rings, which avoids the use of Bessel functions 
and their zeroes. 

It is important to underline that the method that we 
propose is systematic and that it can be applied to rings 
of arbitrary width and curvature. 

The paper is organized as follows: in Section |n| we 
describe the numerical implementation of our method, 
using the Conformal Collocation Method of [1, 2 ; in Sec- 
tion |III| we discuss the extension of the analytical tech- 
niques of [l] to the case of annular domains; in Section 
|IV| we apply the methods, both analytical and numeri- 
cal, to the solvable problem of a circular annulus and to 
the less tractable problem of a "Robnik's ring", i.e. a 
ring whose external border corresponds to the family of 
quantum billiards studied by Robnik in [7] and known as 
"Robnik's billiards". Finally in Section |V| we summarize 
our results and draw our conclusions. 



II. CONFORMAL COLLOCATION METHOD 

In this section we describe the application of the Con- 
formal Collocation Method (CCM) of Ref. [1] l2J to the 
solution of Schrodinger equation in a quantum ring. Al- 
though this approach has been already described in detail 
in those papers, we briefly review it here to make the dis- 
cussion self-contained and to highlight the modifications 
which are needed to implement the specific problem at 
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FIG. 1: (color online) Annulus obtained conformally mapping 
a rectangle of sides 2Lx — 2 and 2Ly — 27i centered in the 
origin. 



Our starting point is the homogeneous Helmholtz 
equation on a "ring-like" domain P, which can be confor- 
mally mapped to an inhomogeneous Helmholtz equation 
on a "simpler" domain which is assumed here to be a 
rectangle. Let w = f{z) where w = u-\-iv and z = x-\-iy 
{{u^v) e V and {x^y) G Q). 

As a result one is left to work with the inhomogeneous 
Helmholtz equation on Vt: 



1 



(1) 



where S e 

An explicit example of conformal map with the desired 
properties is clearly the exponential map 



(2) 



which maps a rectangle of sides 2Lx and 27r centered in 
the origin into an annulus of radiuses e~^^^ and 1 re- 
spectively. Fig. [T] displays the annulus obtained mapping 
a rectangle of sides Lx = 2 and Ly = 2i\ centered in 
the origin (clearly using a smaller Ly one would obtain 
an arc instead of the full ring). Notice that the func- 
tion of eq. ([2| maps the horizontal sides into the hori- 
zontal segment which cuts the ring on the negative axis. 
This is clearly a different situation from those considered 
in Ref. [1 : as a matter of fact in this case the direct 
approach of Ref. [1 would describe an annulus with a 
straight cut on which Dirichlet boundary conditions are 
obeyed. Although this is also an interesting problem in 
itself, we want here to treat a ring with no cuts. 

The solution to this problem is obtained by imposing 
periodic (Dirichlet) boundary conditions on the horizon- 
tal (vertical) sides of the rectangle. Under the conformal 



map of eq. ([2| the rectangle is then mapped into an annu- 
lus fulfilling Dirichlet boundary conditions on the smaller 
and larger circles. To implement the new boundary con- 
ditions in the problem we thus need to introduce a proper 
set of Little Sine Functions (LSF): this set corresponds 
to LSFx of Ref. [8]. 

We report here the explicit form of the LSF\\ 
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with k = — 7V/2, . . . , A/'/2. Notice that we use here a 
different convention for N with respect to Ref. [8 . These 



functions define the homogeneous grid Xk 
satisfy the orthogonality relation 



2Lk 



j si'\N,L,x)sf\N,L,x)dx ■ 



2L 



7V + 1 



and 



(4) 



Notice that h^^"^ = is the grid spacing of the LSFi. 

The LSF2 fulfilling Dirichlet boundary conditions may 
be cast in the form 



si''\N,L,x) = 



(-1)^ cos (f) sin (^) 
N sin (If)- sin (f)' 



(5) 



with A; = -7V/2 + l,..., N/2 - 1 ^ In this case these func- 

2Lk 

tion define the homogeneous grid Xk = and obey the 
orthogonality relation 



J-L ' 



si--> {N, L, x)sf^ {N, L,x)dx = =^5kj • 



2L. 
AT" 



(6) 



^ is the grid spacing of the LSF2. 



Notice that h'^^^^ 

For a given N (even integer) there are + 1 (A^ 
1) LSF functions obeying periodic (Dirichlet) boundary 
conditions, each peaked (with value 1) at a point Xk and 
vanishing at the remaining grid points Xj, j 7^ k. 

A function f{x) obeying periodic be may be interpo- 
lated using the [h^N^x) as 



N/2 



k=-N/2 



(7) 



Similarly we may derive twice this expression to obtain 



k=-N/2 



dx'^ 



^ The reader may check that the expression given here for these 
function is equivalent to the one used before in Ref. jlj 
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s^f\h,N,x) 



k^-N/2 j^-N/2 
N/2 N/2 

E E fixk)cif sf\h,N,x), (8) 

k^-N/2 j^-N/2 



where in the last hne we have introduced the matrix 
which provides a representation for 



the second derivative operator on the grid. The case of 
Dirichlet be has already been discussed in Ref. Hj, al- 
though it can be obtained straightforwardly repeating 
the same steps done here. To take into account the pres- 
ence of different sets of LSF we modify the notation of 
Ref. P and call c^.^^'' the matrix elements of the second 
derivative obtained with Dirichlet be. Notice that the 
latin indices /c, j span different values in the two cases 
and that the grid points also differ in the two cases. 

Omitting some trivial steps (see Ref. [l] for more de- 
tail) we may now may easily discretize eqn. ([T]) using 
Dirichlet be in the x direction and periodic be in the 
y direction. A suitable "basis" for this discretization is 
obtained with the direct product of the LSF in each di- 
rection: 

^ ^sf\h,N,x)s'^i){h,N,y) = -Y, ^ 



ir^-^fc'i' sf'\h,N,x)s'^l\h,N,y) (9) 

where it is understood that the latin indices span different 
ranges for the two LSF. 

To obtain the matrix element of the operator on the 
grid we need to associate a single integer to any pair of 
indices which define an element of the two dimensional 
grid. We may write: 



k 
k' 



K-{N-1) 
N 



K 



7V-l + e 



TV 



K 



N-l 



(10) 
(11) 



where e ^ 0^ and [a] means integer part of a. In this 
way we are able to identify a point of the grid in terms of 
a single integer which takes values from 1 to A^^ — 1. 

Using these relations we may read off the matrix ele- 
ment of O as 



O 



KK' 



Sim 

-kj 



SkjC 



{12) 
k'j' 



(12) 



Although the procedure described above uses grid with 
the same N for the x and y direction, a more appropriate 
choice is to use meshes with the same grid size on each 
orthogonal direction, i.e. h^^^ ~ h^^^\ In this way we 
may establish the relation between the number of grid 
points on each direction 



AT, 



{Ny + 1) . 



We may easily understand the physics contained in this 
relation: for thin rings, Lx <C Ly, the number of grid 
points in the x-direction is much less than the number 
of grid points in the ^/-direction, since the excitation of 
trasverse modes requires much higher energy than the 
excitation of longitudinal modes. 

In this case one can extend the previous relations for 
the grid to: 



k = K-{N,-1) 

k' - 



K 



A". 
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K 



A^.-l + e, 



The implementation of these considerations allows us 
to represent the differential operator as a (A^^^ — '\-){Ny + 
1) X [Nx — ^){Ny + 1) hermitean matrix. 

It is useful to summarize the differences with the re- 
sults of Ref. [1 : 

• The region 1] is a rectangle of sides 2Lx and 2L^, 
where L^y = tt to allow a closed ring 

• The operator O is represented on the grid by a 
(A^^ - 1)(A^^ + 1) X {Nx - l){Ny + 1) hermitean 
matrix 

• The collocation points, corresponding to the nodes 
of the LSF functions, differ for the x and y direc- 
tions: on the X direction they are distributed fol- 
lowing the zeroes of the LSF with Dirichlet bound- 
ary conditions, while on the y direction they are 
distributed following the zeroes of the LSF with 
periodic boundary conditions. Their numbers (A^^^ 
and Ny) are also different, although the grid size is 
(approximately) the same. 

Apart from these small differences, the remaining fea- 
tures of the collocation method are unchanged. In par- 
ticular, the matrix representing O on the grid is obtained 
from the product of a diagonal matrix, representing 1/E 
on the grid, with a non-diagonal sparse matrix, repre- 
senting the Laplacian with mixed be on the grid. While 
the first matrix is specific to the problem considered and 
therefore it needs to be calculated each time that a differ- 
ent shape is chosen, the second matrix is universal and 
therefore it can be calculated and stored once and for 
all. As before the calculation of the diagonal matrix is 
not computationally demanding since it just requires the 
evaluation of the function 1/E at the {Nx — l)(A'?y + 1) 
points forming the grid. 

We will illustrate later specific applications of the 
CCM. 



III. ANALYTICAL METHODS 

In this section we generalize the analytical approach of 
Ref. [1] to describe quantum rings. We may consider two 
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different approaches: in the first one the quantum ring 
is obtained performing a conformal map of a rectangle 
centered in the origin (as described in the previous sec- 
tion); in the second one the quantum ring is obtained by 
applying a conformal map directly to a circular annulus. 

While the formulas of Ref. [1 hold both cases, one 
needs to work with different orthonormal basis in the 
two cases. 

In the first case where 1] is a rectangle of sides 2Lx 
and 27r, the basis is obtained by the direct product of 
functions obeying Dirichlet and periodic boundary con- 
ditions, i.e. 



1 



. ( , T \ 



2L, 



for the Dirichlet be and 

Xo(^) = 

Xny{y) = 



27T 



: COS (uyy) , 



(13) 

(14) 
(15) 
(16) 



(l^riyiy) = sin (uyy), 

for the periodic be. Notice that rix^y = 1,2, — 
The basis on Q may then be written as 



where the value Uy = can only be reached for the states 
with 5 = 1. 

In the second case, the wave function of a circular an- 
nulus with a < r < b is (see for example Ref. 0) 



{r,e)=Nr. 



COS nO 
sin nO 



a 

s = 1 
s = 2 



where Jm and Ym are Bessel functions of first and second 
kind and k is the n^^ root of the equation 



YjYi (fc) JfT 



Jm (^) 



kb 



a 







(17) 



The energies of the annulus are then given by 

2 



a 



, m = 0,1,2,... , n = 1,2... (18) 



Although one may work equally well with each of the 
two basis, from the point of view of an analytical calcula- 
tion the first one offers the advantage of simplicity, since 
it involves only elementary functions. We will therefore 
focus on this basis. 

We consider the most general conformal transforma- 
tion which maps the rectangle Vt onto a ring of arbitrary 
shape: 



C -g{z) , 



(19) 



k=0 



where r^o = 1 and C > is a constant factor, represent- 
ing a dilation. The energies of the ring obtained using 
the mapping g{z) (En) are related to those of the ring 
obtained using g{z) {En) by the simple relation ^ 



En. = 



En 
C2 • 



We may therefore work with g{z) and then simple rescale 
the energies obtained. 

Notice also that with the choice C = 1 and r]k>o = 
we obtain the map to the circular annulus considered 
earlier. 

With simple algebra we obtain the conformal density 



oo oo 



fe=0 j=0 
N(fe+i)+2 



^^.-L.).-..- cos(y(fc-j)) 



(20) 



and cr(x, y) = i;(x, y) — 1. 

As discussed in Ref. [1 we need to calculate the matrix 
elements of a between the states of We find: 



{nx,ny,s\a{x,y)\n'^,ny,s') = -6n^^n'Jny,n'Jss' 

oo oo 

"^"^VkVjik + l)(j + l)An,^ny,s,n'^Xy,s'{kJ^21) 

k=o i=o 



where 



nx,ny,s,n'^,n'y,s 



/-\-Lx P + 

dx I dy ^nx^nyA^^y) 
-Lrr J —TV 



COS (y(fc-j)) (^'V) 



= Q;„„ a„; Ss.s' W„„ (fc, J )4'^,„; ,k,. 



(22) 
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and Qfo = 1/\/2ti5s^i and an^>o = l/v^- 



We have also introduced the definitions 



and 



ny,n'k,j 



r(2) 

ny,n'k,j 



X 

8n^L^n^n',{j + k + 2)e-2^-(^+'=+2) (^(_i)»x+<e2i-(^+'=+2) - 1 
(4L2 (i + A; + 2)2 + 7r2(n, - (4L2(j + + 2)2 + tt 2(n, + 



/+7r 
cos{nyy) cos(n'yy) cos{{k - j)y)dy 
-7T 



/+7r 
sin(n^?/) sin(n^?/) cos((A: - j)y)dy 
-TT 



^n^-n^ + (/c-j) ^ny-\-n'y-\-{k-j) ^ny-n'y-{k-j) ^ny+n'y-{k-j) 



(23) 



(24) 



Notice that the t[^\, ju have been already defined in are the "unperturbed" energies (i.e. the energy on the 
Ref. [Ij; the reader should also observe that there are no rectangle 
terms mixing states with different values of s and that 
4'in;,fc,, = ^nyXy,k,j ™less + n; = ±{k - j). 



This expression has been derived in Ref. \V resumming 
We may understand the physical consequences of these ^P^^^^*^ ^erms in the perturbative expansion 'correspond- 



properties by considering the approximate expression for l^S. geometric series: it has later been applied in 

Ref. \iO\ to obtain a Weyl-like law for inhomogeneous 
drums. 



the energy: 



where 



s\'E{x,y)\n^,ny,s) 



(25) 



(26) 



With simple algebra we obtain 



00 00 



{n^,ny,s\^{x,y)\n^,ny,s) = ^^r]kVj{k ^ ^ 1) [l^n,,n,,fc,fc4j + (-l)^^^^n,,n,,/c,j(^2n^-|/c-j|] 

k^Oj=0 



J2vl{k + lfWn^,n.,k,k + (-1)'+' Yl VkVk-2ny{k + l){k + 1 - 2ny)Wn^ ,k,k-2ny 
k—0 k—2ny 
oo 

,71^ ,k,k-\-2ny 1 

(27) 



k=0 



where the terms depending explicitly on s will be respon- 



sible of the breaking of the degeneration of levels. 
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IV. APPLICATIONS 



We will now consider few applications of the numerical 
and analytical techniques developed here and in Ref. [I] 
to quantum rings of different shape. 



A. Circular annulus 

Our first application is to the calculation of the ener- 
gies and wave functions of the circular annulus using the 
basis of the rectangle. 

We may use eq. (25) and obtain an explicit formula for 



the energies of a circular annulus: 

2 (log^(a) + TT^n^^) {uy^ \og\a) + tt^ 



TT^ (a^ — 1) Ux'^ log(a) 



where a = e~^^^ is the inner radius of the ring. 

Notice that the states corresponding to = are non- 
degenerate, whereas the states corresponding to n^y > 
are doubly degenerate: in other words the basis that we 
are using reproduces the exact pattern of degeneration of 
the circular annulus. 

In Fig. [2] we compare the exact numerical results ob- 
tained solving eq. ( 18 ) (solid curve) with the results ob- 



tained using the analytical formula (28) (dashed curve) 
for the first 2000 states. The dotted line corresponds 
to the numerical results obtained using the CCM with a 
grid with = 14 and Ny = 400 (remember that the 
ratio Nx/{Ny + 1) should approximately be Lx/Ly). Fi- 
nally the dot-dashed line are the results obtained using 
Weyls'law supplemented by Weyl's conjecture 



En 



47rn 



L Uirn 



(29) 



where A and L are area and perimeter of the ring respec- 
tively (Dirichlet boundary conditions are chosen). 

Amazingly the first three curves are practically indis- 
tinguishable. In Fig. [3] we display the the error of the 
analytical formula, which is about 0.1 % for all the states 
considered. 

We may also obtain precise results for the ground 
state of the ring using a variational approach. Follow- 
ing Ref. [1 we pick a trial state 



(30) 



where \k) = \kx,ky) are eigenstates of the rectangle with 
mixed be. Using the inverse operator we may gen- 
erate a new state 



l*^'^) = o-i|x) 



(31) 



which has a larger overlap with the true ground state of 
the system. 
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FIG. 2: (color online) Energy of the first 2000 states of an 
annulus with a — 9/10 and 6=1. The solid line are the exact 
results of eq. |T8| ; the dashed line corresponds to eq. (28); 
the dotted line are the numerical results obtained with CCM; 
the dot-dashed line are the results obtained with Weyl's law 
supplemented by Weyl's conjecture. 
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FIG. 3: (color online) Relative error of the analytical formula 
( 28 ) with respect to the exact results 



The expectation value of O in this state is 



E, 



(32) 



where 



^ ^(0)^(0) 



k,l,'m=0 



7 



.(0)J0) 



= V 

l,l,m,m—0 

X {^|Ei/2|m) . 



-{m'|Si/2|O(«'|S|0 
(34) 



By minimizing Eq^"^ with respect to the coefficients c^^^ 
one can now find precise estimates for the ground state 
energy. Let us now pick a specific trial state. We wih 
work in the hmit of a annulus with {b — a) <C 1, which 
corresponds to using <^1. 

We use the trial state 

Ix) = 4°Vi(^)Xo(^) 

N 



neglecting in first approximation the excited states in the 
X direction. 

We need to calculate the matrix elements 



(l,m,|Ei/2|l,g 
/ " ^/;?(a;)e^-'^-da; 

J Lnr 



■ sinh(L^) 



(m|E|0 



{l,my\i:\l,ly) 



(36) 



(1 



-4Lx 



4(4L2+7r2) 



(37) 



Therefore 



(1) 



' sinh(La; 



1 2 



^2q-Lx sinh(La;) 



7r2(l-e-4^ 
4(4L2+7r2 



Therefore we may write the matrix elements: 

{m\Y}'^\l) = (m,,0|Ei/2|z^,o) = 



7r^(l- 



(38) 



where 



TT 

4L2 



(39) 



where [|] is the integer part of k/2. Since one of the 
coefficients c^^^ can be chosen arbitrarily we may pick 
^^^^ — 1 and minimize E^^ with respect to the remaining 

, for /c > 1, which is 



(0) 



coefficients. The result is c 
equivalent to working with eq. (25). This fact should 
not be surprising, since the states corresponding to non- 
vanishing c^^^ (for k >1) break the rotational invariance 
of the true ground state of the annulus and therefore 
they must raise the energy. This explains the very good 
precision of eq. (28) for annuli with b — a ^ 1. 

To improve this formula for states with smaller inner 
radius we may form a variational trial state with a su- 
perposition of wave functions in the x direction: 



AT 



Ix) = My) E 



(40) 



/. 



^''m'^ sinh(La;) 



x — Lsi 



dx 



Ll+TT^Lxml 
^^L^m,,^,,(e2^^-(-l)'^^+^^) cos(7r(m^+^^)) 

16Lj+87r2L2(m2+^2)+^4(^2_^2)2 

r+Lx 



^mx 



dx 



Lx 



167r^e" 



4(4L3+7r2L^m2) 
-I^.^xm.(e^^^-(-l)^^+-^ 



256Lj+327r2L2(^^2^^^2)+^4(, 



1 ^x ^a; 
I cos{7r (I x+rux)) , , 

1 "^X ^X 



(41) 



Notice that in this case 



4L2 • 
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In Table [T] we report the ground state energy of the 
annulus calculated using the variational principle, i.e us- 
ing the variational formula (38) together with eq. (40), 



with different numbers of variational parameters. The 
case = 1 corresponds to the simple analytical formula 
( [28] ), which is seen to work very well for larger values of 
a. It is easy to understand why more and more terms are 
needed as a gets smaller and smaller: since L gets larger, 
the conformal map strongly deforms the radial coordinate 
in the annulus. Looking at Fig. [l] which corresponds to 



a = we observe that the radial grid spacing is finer 
for smaller values of r. Therefore one needs more terms 
in eq. ([38| to obtain good estimates of the energy. Notice 
that this argument also holds for the Conformal Collo- 
cation Method, since also the collocation points are also 
distributed non-uniformly in the radial direction, as the 
central hole is made smaller. 



a 


N = 1 


N = 2 


N = 3 


N = 4 


N = 5 


N = 10 


exact 


1/100000 


24.74037682 


13.58247813 


10.13874976 


8.585413233 


7.763566890 


6.639748109 


6.493374419 


1/10000 


20.56383831 


11.77986914 


9.129597650 


7.994645335 


7.425906261 


6.751169408 


6.693191961 


1/1000 


16.67307038 


10.21973107 


8.409275276 


7.691356864 


7.367206060 


7.062604524 


7.048038125 


1/100 


13.49800490 


9.359102734 


8.360288280 


8.036684452 


7.920642510 


7.846777753 


7.845162680 


1/10 


13.31090873 


11.31798812 


11.04543198 


10.99627372 


10.98600010 


10.98222952 


10.98218981 


1/2 


39.81860401 


39.03110135 


39.01558313 


39.01353082 


39.01336954 


39.01328910 


39.01328850 


9/10 


987.1575032 


986.6840018 


986.6837734 


986.6831548 


986.6831530 


986.6831306 


986.6831305 


99/100 


98696.22334 


98695.79206 


98695.79206 


98695.79151 


98695.79151 


98695.79149 


98695.79149 



TABLE I: Energy of t he ground state of the annulus with varying inner radius calculated using the variational formula (38) 
together with eq. (40), with different numbers of variational parameters {N — 1 is the number of variational parameters) . 



Notice that in this case 



In Fig. [4] we have plotted the exact energy of the 
ground state of the annulus as a function of the inner 
radius a and we have compared it with the different ap- 
proximations considered earlier. The simple analytical 
formula of eq. (28) is seen to work quite well for thin 



annuli, while the most accurate variational calculation of 
Table [l] is seen to reproduce the exact results even for 
very small inner radiuses. We also report the CCM re- 
sults corresponding to a grid with = 50. 



B. Robnik's rings 

We consider the conformal map 



(42) 



which maps the rectangle of sides 2Lx and 27r into a 
deformed ring, which for a = has the shape of a circular 
ring and for a = 1/2 has the shape of a cardiod ring. We 
have called the family of these annular billiards "Robnik's 
rings", in analogy with the simply connected billiards 
known as "Robnik's billiards" having the same external 
contour [7 . 

In this case we have: 



4ae 



3ic 3-£/^ 



cos{y) 



(43) 



Notice that in this case the width of the ring is not 
uniform: we call the smallest and largest width of the 
ring respectively, which read 



r± = ±a(l-e"'^^-) -e"^^- + 1 

Their average is insensitive to a: 

_ _ r+ + r_ 



1 - e 



(44) 



(45) 



Notice also that 1 — f = e~^^^ = a (in this case a is the 
average inner radius of the quantum ring). 

Using the general results obtained earlier we may write 
the analytical formula: 



2 (log^(a) + TT^n^^) (n^^ log^(a) + tt^h^^ 
7r2 (a2 _ 1) log(a) + a^R„^ (a) 



where 



RnA^) = 2^2 (a^ - 1) log(a) (log'(a) + ir^n,^) . 



This formula reduces to eq. ([28]) for a = 0. 

In Figj6]we display the energies of the first 2000 states 
of an annulus obtained conformally mapping a rectangle 
of sides 2Lx = 1/5 and 2Ly = 27r with the map f{z) = 
z + az'^ using a = 1/10. The solid line represents the 
precise numerical values obtained using CCM with a grid 
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FIG. 4: (color online) Ground state energy of the annulus as 
a function of the inner radius keeping the outer radius fixed 
(6 = 1). The dashed line is the exact result, while the dashed 
line is the analytical formula of eq. (28). The dotted line is 

the pluses are 
10. 



the result obtained using CCM with = 50; 
the variational results of Table [l] corresponding to N 



FIG. 6: (color online) Energies of the first 2000 states of an 
annulus obtained conformally mapping a rectangle of sides 
2Lx = 1/5 and 2Ly = 27r with the map f(z) — z ^ az^ using 
a = 1/10. The solild line are the numerical results obtained 
using CCM; the dashed line corresponds to the analytical for- 
mula, the dotted line is Weyl's law. 




FIG. 5: (color online) Asymmetric annulus obtained confor- 



mally mapping a rectangle of sides 2Lx 
centered in the origin, a = 1/10. 



1/5 and 2Ly = 27i 



with Nx = 14 and A^^ = 400; the dashed line corresponds 
to the results of the analytical formula of eq. (46); finally, 



the dotted line corresponds to Weyl's formula eq. (29). 
In Fig. [6] we display a detail of the previous figure, for 
highly excited states around n = 2000. 

Looking at the figure we see that our analytical formula 
describes more accurately than Weyl's equation, eq. (29), 



the low lying part of the spectrum, although the latter 
describes better the high part of the spectrum. We also 
notice that the analytical formula displays a behaviour 
already observed for the circular annulus: around spe- 
cific n, corresponding to the opening of a new trasversal 
mode, the curve displays a small change of slope. While 
this behaviour is correct for the circular annulus, where 
it is indeed observed, neither the numerical results nor 
(of course) Weyl's expression display such a behaviour. 
However we may easily understand the origin of it: the 
analytical formula eq. (46) is obtained using a resum- 



mation of the terms in the Rayleigh-Schrodinger series 
which contains only the expectation values of the confor- 
mal density in a given state. For the circular annulus, this 
approximation is quite good because there is no mixing 
between transversal and longitudinal modes as a result 
of the rotational symmetry of the domain: this is re- 
flected in the amazing accuracy of the analytical formula 
in this case. For a general annulus, the domain is not 
invariant under rotation and transverse and longitudinal 
modes actually mix. Already to second order in pertur- 
bation theory, there are terms which allow this mixing, 
although these are not present in the resummed formula. 

As a check of the quality of the numerical results ob- 
tained with the CCM we have applied a method devised 
by Berry, ref. [11], which allows one to extract the geo- 
metrical feature of a domain from a limited sequence of 
eigenvalues, using improved eigenvalues sums. For the 
area, for example, one obtains a family of approximating 
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FIG. 7: (color online) Detail of the previous figure. 



functions whose first few elements read [TT] 



N 



n=0 
N 

Ai{t) = 47rt'^e-^-\2Ent-l) 

n=0 
N 

A2{t) = 47Tt^e-^^^Ent{2Ent - 3) 



(47) 
(48) 
(49) 



n=0 
N 



A3{t) = 47rt ^ e-'^"*^„t(4£;2t2 _ i2^„t + 3)/^50) 

(51) 



n=0 



Analogous expressions are given in Ref. [TT] for the 
length and constant terms in the asymptotic expasion of 
the function counting the number of states up to a given 
energy [12] 



{n{E)) ^ AE/47T - lVe/47t + C ■ 



(52) 



The constant C should vanish for domains with one hole. 

In figs |8|9| and we have displayed the approximants 
for the area, perimeter and for the constant term built 
using the first 2000 numerical eigenvalues obtained with 
the COM with a grid 14 x 400. The horizontal lines are 
the exact values. 

To estimate the optimal value of t where the ap- 
proximant should be evaluated we minimize the sum 
of the squares of the derivatives of each approximant. 
For example, for the area we minimize the quantity 
A[{t)'^ ^ A2{t)'^ ^ A'^{t)'^ (we leave out Ao{t) because it is 
clearly less precise). For the length and constant terms 
we use the approximants displayed in the figures. 

As shown in Table [llj the method of Berry allows us 
to extract the area, perimeter and constant terms of the 





A 


L 


C 


method of Ref. 1111 
exact 


1.07027 
1.07032 


11.5221 
11.5250 


-0.01738 




TABLE II: Estimates for the geometrical features of the 
annulus obtained conformally mapping a rectangle of sides 
2Lx = 1/5 and 2Ly = 27v with the map f{z) = z-\- otz^ using 
a = 1/10 using the method of Berry. [TT] on the numerical re- 
sults obtained with the COM. The last row contains the exact 
results. 




0.01 



FIG. 8: (color online) Area approximant functions Amify for 
the first 2000 eigenvalues obtained with CCM with a grid 
14 X 400. 



annulus with a truly amazing precision. Notice that the 
exact values of the area and perimeter are obtained di- 
rectly using the conformal density: 



A = dxdy i;(x, y) 
Jn 

L = ds ^/^{x,y) , 
Jan 



(53) 
(54) 



where Q is the rectangle of sides 2Lx = 1/5 and 2Ly = 2i\ 
and dVt is the border of Vt where Dirichlet be are obeyed ^ . 

As we did before, we may also resort to a variational 
calculation and use the trial state: 

Ix) = 4°Vi(^)Xo(y) 

N 

+ ^i(^) E Ki-i'^n.(y) + 4°U»«M (55) 



neglecting in first approximation the excited states in the 
X direction. In the case of the circular ring we have seen 



^ In ref. 10 these relations are used to derive Weyl's law from 
perturbation theory. 
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FIG. 9: (color online) Length approximant functions Lmif) 
for the first 2000 eigenvalues obtained with COM with a grid 
14 X 400. 
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FIG. 11: (color online) Energy of the ground state of the 
asymmetric annulus with Lx — 1/10 and a — 1/10. The 
horizontal line is the precise numerical value obtained with 
the CCM with a grid with 14 x 400; the pluses are the vari- 
ational results obtained with a varying number of variational 
parameters. 
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FIG. 10: (color online) Constant approximant functions 
Cm{t) for the first 2000 eigenvalues obtained with CCM with 
a grid 14 x 400. 



that only the first term contributes, due to the rotational 
invariance of the ground state; this limit corresponds to 
the simple formula (25), which was seen to work remark- 



ably well for thin circular rings. 

In the present case we do not expect that the lowest 
order formula work well even for thin rings, due to the 
varying width of the ring, which favors the breaking of 
the rotational invariance. 

We may apply this variational calculation to the Rob- 
nik annulus that we have considered before, which corre- 
sponds to a = 1/10 and = 1/10, shown in Fig. [sj In 



this case we have 



= 0.148 



0.214 



(56) 



In Fig. 11 we show the variational estimates for the 
ground state energy of this ring obtained using a vary- 
ing number of variational parameters Nyar (the case 
N„ar = corresponds to the simple analytical formula 



(25)); the horizontal line is the precise value obtained 
with the Conformal Collocation method with a grid with 
N = 100. As anticipated, the eq. (25) in this case pro- 
vides a poor approximation to the exact energy, lead- 
ing to a substantial overestimation of it. On the other 
hand, the inclusion of few longitudinal modes leads to 
quite drastic improvement, which can be already appre- 
ciated for Nyar = 2. We should point out that, even 
letting Nyar => oo (assuming that we can perform the 
calculation), we do not expect that the variational en- 
ergy converge to the exact energy: the reason of this is 
that the variational ansatz that we are using includes 
only longitudinal modes (in the y direction) while com- 
pletely neglecting the trasversal ones (in the x direction) . 
For thin rings, the latter are less important, and there- 
fore we are able to obtain good (but not arbitrarily good) 
approximations; for thick rings, on the other hand, only 
the inclusion of both modes will lead to acceptable ap- 
proximations. 

Finally in Fig. 12 we have plotted the wave function 
of the ground state of this annulus: as anticipated, this 
wave function is not invariant under arbitrary rotations 
and it is peaked in the region where the ring is wider, as 
expected. The wave function is multiplied by a factor 10 
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FIG. 12: (color online) Wave function (multiplied by a factor 
of 10) of the ground state of the Robnik's annulus with Lx = 
1/10 and a = 1/10. 



to better appreciate the localization in the region where 
the width of the annulus is larger. 

V. CONCLUSIONS 

In this paper we have extended the powerful techniques 
of Ref. [1, 10 to the study of quantum rings; working 
on specific examples we have obtained both numerical 
and analytical results of remarkable precision. In fact 
our approach provides a systematic tool for the study of 
the spectrum of quantum rings of arbitrary shape, not 
relying on specific assumptions concerning the shape or 



the width of the ring. Clearly, the problem of finding 
the particular conformal map which sends the rectangle 
to a specific ring may be attacked numerically, when the 
explicit expression is unknown [I3l [14] . 

The results that we obtain in the paper maybe roughly 
classified into: 

• numerical results obtained using the Conformal 
Collocation Method (CCM) for a large number of 
levels {'^ 2000); the accuracy of these results has 
been verified using Berry's approximants to esti- 
mate the geometrical features of the rings with a 
selected number of levels, calculated numerically; 

• variational estimates for the ground state of a quan- 
tum ring; 

• explicit analytical formulas for the whole spectrum 
of a quantum ring; the relation of these formulas 
with Weyl's law is also obtained in light of the re- 
sults of Ref. [10 ; 

We believe that the importance of our methods thus 
relies in the possibility of describing with precision both 
the low and high end of the spectrum of a quantum ring 
and on the possibility of improving this precision in a 
simple and sistematic fashion. To the best of our knowl- 
edge there are no other approaches in the literature which 
share these features. 

Finally, we wish to mention that the approach that 
we have described here can be applied directly also to 
annular drums of variable density (for simply connected 
drums this has been done in Ref. |10j). 
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